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ZmAMVi 

In  this  pctper  tne  funct Iona l-e..|Uat ion  approach  of  aynamic 
progi'aTjmlng;  is  used  to  trf-at  a  iineai  programming  problem 
lnvolvlr\e  a  elightly  IntertAlneJ'  matrix — l.e.,  one  that  is 
almosr  hiocK  dlagci»al . 


1 .  INTRODUCTION 

Consider  the  problem  of  maximizing  the  linear  form 
3N 

(1)  L.(x)  .  Z  X 

^  1-1  ^ 

over  all  satisfying  the  constraints 
* 

(2)  +  &i2^2  ^  ‘^l* 

^21^1  ^22^2  ^  ^25^3  ^  ^2* 

®3l^l  ^33^3  ■*“  ^1*4  -  ^3' 

+  ^43^5  ^  ^46*6  ^  ^4* 

^54^4  +  +  ^56^6  - 

^64^4  ^  ^63^5  •  We  ^  ^2^7  "6' 


^3N-2,3N-2^3N-2  ^3N-2,3N-1^3N-1  ^  ®3N-2,3N^3N  °3N-2' 
^3N-l,3N-2^3N-2  ^3N-1,3N-1*3N~1  ®3N-1,3N*3N  ^  ^3N-1' 
®3N,3N-2^3N-2  ®3N,3N-1^3N-1  ^3N,3N^3N  ^  ^3N'  . 

and 

(3)  2  0- 

It  Is  assumed  throughout  that  a^j  2  *^1  ^  ®1  ^  with 

a  sufficient  set  of  the  a^j  positive  so  that  the  maximum  of 
Ijj(x)  Is  not  Infinite. 


We  wish  to  attack  this  problem — which  arises  from  the  study 
of  weakly  coupled  economic  systems,  or  alternatively  from  the 
study  of  multistage  processes  with  almost— Independent  stages— 
by  means  of  the  techniques  of  dynamic  programming.  Specifically, 
we  shall  show  that  the  computational  solution  can  be  obtained 
by  means  of  a  sequence  of  functions  of  one  variable. 


It  will  be  clear  that  a  number  of  similar  problems  Involving 
0 

almost  bl-ock-dl agonal  matrices  can  be  treated  by  means  of  the 
same  general  method.  In  another  paper,  f 2 J ,  we  have  Illustrated 
the  application  of  this  Idea  to  mechanical  and  electrical  systems 
In  which  the  matrices  are  symmetric. 


2.  DYNAMIC  PROGRAMMING  FORMULATION 

Let  us  define  the  sequence  of  functions  of  z, 

(1)  fQ(z)  -  0,  f^(z)  -  Max  ly.{x),  K*1,2,...,N, 

where  the  are  subject  to  the  constraints  given  above,  with 
the  exception  that  the  last  constraint  Is  now 


®3K,.^K-2^5K~2  ^:5K,3K-1*:5K~1  ^  ^ 


Employing  the  principle  of  optimality,  cf.  [l]  ,  we  see  that 
the  sequence  satisfies  the  recurrence  relation 


Max 

[^3K-2'*5K-1'^5k] 


^3K-2  '*■  *5K-1  ^K-i 


"  ^K-1^5K-2^  J' 

with  the  variables  constraints 


^5K-2,5K~2^3K-2  ®3K-2,3K-1^5K-1  ®3K-'2,^K^5K  ^  ®3K-2' 

®3K-1»3K~2^3K-2  ^  ®3K-1,3K-1^5K-1  "*■  *5K-1,3K*3K  ^  ®3K-1' 
®3K,3K-2^3K-2  ^  ®3K,3K-1*3K~1  ®3K,3K*3K  ^ 

^K-1^3K-2  ^  '=3K-3' 

^3K-2 '^3K-1 '’'3K  ^ 

3»  SIMPLIFICATION 


We  can  write 

(1)  Tj^iz)  « 


where  Rj^  Is  the 


the  recurrence  relation  of  (2.3)  In  the  form 


Max 

Max 

[• 

..]  ] 

^3K-2 

^  *3K-1  '*3K 

Max 

Y 

Max  (x 
L  R 

?K-1 

X^j^)  +  X^j_, 

*3K-2 

^  ‘■k-1 

^®3K-J 

^  ^K-1^3K-2 

region 

In  (J'}k_i>>'3k) 

space 

defined  by 

®3K-2,3K-1^3K-1  ®3K-2,3K^3K  ^  °3K-2  “  ®3K-2,3K-2^3K-2' 
^3K-1,:JK-1*3X-1  ^  ®3K-1,3K^5K  ^  ^3K-l  ~  ®3K-1 ,3K-2^3K-2 ' 
®3K,3K-1^3K-1  ^  *3K,3K*3K  ^  ~  ^3K,3K-2^3K-2 ' 


(2) 


